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ANOMALOUS  TRANSPORT  FROM  PLASMA  WAVES 


I 


Can  one  derive,  from  first  principles,  simple  selling  laws  for  complicated  plasma  devices 
like  a  lokamak'’  For  instance  the  current  vogue  is  to  scale  confinement  in  (okamaks  linearly 
with  density,  so  clearly  one  would  like  a  thermal  conduction  coefficient  scaling  as  n  The 
problem  is  that  this  thermal  conduction  must  he  derived  from  the  quasi-linear  and  nonlinear 
theory  of  very  esoteric  instabilities  tven  if  someone  could  take  these  instabilities  and  derive  a 
transport  coefficient  which  t'jd  all  the  right  magnitudes  and  scalings,  would  anyone  else  either 
understand  it  or  believe  it'* 

My  own  opinion  is  that  one  cannot  derive  simple  scaling  laws  in  an  understandable, 
believable  way  How  then  docs  one  explain  or  predict  the  behavior  of  complicated  plasma  dev¬ 
ices  which  arc  dominated  by  anomalous  transport '  One  possibility  of  course  is  to  use 
phenomenological  transport  coefficients,  and  this  may  in  fact  work  well  It  is  certainly  likely 
that  current  devices  could  be  extrapolated  a  factor  of  two  in  every  parameter  this  way  Howev¬ 
er  this  is  somewhat  unsatisfying,  one  would  somehow  like  to  relate  performance  to  fundamen¬ 
tal.  processes,  if  only  for  fundamental  scientific  reasons  That  will  be  the  focus  of  the 
remainder  of  this  paper 


Let  us  look  at  what  is  needed  to  describe  a  plasma  whose  transport  is  dominated  by  insta¬ 


ll 


bility  As  we  will  see.  there  is  what  might  be  called  a  heirarchy  of  necessary  information  start¬ 
ing  with  the  simplest  questions  of  stability  threshold  and  proceeding  to  the  most  complicated 
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nonlinear  theory  How  high  one  has  to  proceed  in  this  heirarchy  is  then  an  indication  of  how 
easy  and/or  how  reliable  the  theory  is 

Take  for  instance  a  temperature  profile  in  a  mkamak  plasma. and  look  at  some  radius  r„ 
The  most  fundamental  question  is  whether  or  not  ihc  plasma  is  unstable  at  this  point  If  the 
plasma  is  stable  here,  clearly  one  uses  classical  transport  If  n  is  unstable  anomalous  transport 
should  be  involved  Thus  the  first  necessary  piece  of  information  is  the  stability  ihrcshold  As 
we  will  sec.  there  arc  some  circumstances  where  this  is  all  that  is  needed  to  describe  the  sys- 
tern  Since  this  description  only  utilizes  the  first  step  in  the  aforementioned  heirarchy,  it  will 
surely  be  the  simplest  and  most  reliable  theory  of  anomalous  transport  This  is  the  marginal 
stability  theory  which  we  1-4  and  others  5-7  have  previously  discussed 

The  idea  behind  the  marginal  stability  theory  is  the  following  Imagine  a  plasma  which  is 
stable,  but  which  is  forced  by  some  external  mechanism  to  an  unstable  state  As  a  concrete  ex¬ 
ample  let  us  consider  an  ohmically  heated  tokamak  The  current  heals  the  central  region  and 
since  the  edge  is  cool,  an  electron  temperature  gradient  is  forced  upon  the  system  Let  us  now 
postulate  that  when  the  electron  temperature  gradient  exceeds  some  critical  value,  the  plasma 
becomes  unstable  Then  at  this  point,  an  anomalously  large  thermal  conduction  is  generated  by 
the  instability  A  possible  functional  dependence  of  K,.  on  temperature  gradient  is  shown  in 
Fig  1  This  large  Kf  will  cool,  the  plasma,  the  instability  will  shut  off  and  K,  will  go  to  its  clas¬ 
sical  value  However  once  the  plasma  becomes  stable,  the  ohmic  heating  will  increase  its  tem¬ 
perature  and  it  will  again  be  driven  unstable  Clearly  then  there  is  a  dynamic  balance  between 
heating  and  anomalous  thermal  conduction  with  the  plasma  sitting  at  (or  perhaps  oscillating 
about)  the  marginal  stability  point  The  fundamental  quantity  to  determine  then  is  not  the  tur¬ 
bulent  spectrum,  but  the  plasma  profile  Once  we  have  the  profile,  we  can  calculate  the  thermal 
conduction  because  we  know  the  profile  and  input  power  Once  we  have  K,.  simple  quasi- 
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linear  theory  gives  us  the  lluitujtion  lex  cl  Thus  the  logic  of  a  marginal  stability  calculation  is 
just  the  reverse  of  a  conventional  nonlinear  calculation,  as  shown  in  fig  2 


Now  examine  the  second  complication  m  our  heiiarchy  Let  us  say  that  the  marginal  sta¬ 
bility  hypothesis  is  valul,  hut  there  is  not  one  anomalous  transport  effect  but  many  To  contin¬ 
ue  with  our  example  of  the  tokanuk.  n  is  known  that  both  electron  temperature  gradient,  and 
density  gradient  contribute  to  instability  ,  and  that  not  only  is  there  anomalous  thermal  conduc¬ 
tion.  but  also  anomalous  diffusion  In  this  case  marginal  stability  does  not  imply  only  particular 
temperature  gradient  but  rather  some  relation  between  temperature  and  density  gradient  Then 

the  stability  threshold  is  not  sufficient  information  to  determine  the  profile  One  must  also 
know  the  relation  between  k  and  l)  If  one  knows  the  spectrum,  and  if  quasi-lincar  theory  is 
valid,  then  one  can  determine  this  relation  However  now  more  information  is  needed  than  just 
stability  threshold,  correspondingly  the  theory  is  on  somewhat  shakier  ground  than  it  is  if  there 
is  only  a  single  anomalous  effect 

1  et  us  illustrate  how  a  marginal  stability  works  in  a  one  or  two  parameter  space  F  irst  im¬ 
agine  that  only  temperature  gradient  drives  instability  and  thermal  conduction  is  the  only 
anomalous  process  Then  the  steady  stale  equations  arc  heat  balance 


and  Marginal  stability. 


n 


b  ..  9T 

dx  dx 


♦  S; 


(it 


o  --  -  P  -  0  (2) 

dx 

dT 

where  S.  denote  the  energy  sources  In  Eg  <21.  ot-  i*  the  growth  rate  which  is  driven  by 

temperature  gradient  and  P  denotes  the  damping  from  all  sources,  for  instance  shear,  ion-ion 
collisions,  etc  The  marginal  stability  approach  then  consists  of  using  tq  (2)  instead  of  Eq  (1) 
as  an  equation  for  temperature  Equation  (II  is  used  instead  as  an  equation  for  K.  and  it 
reduces  to 
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d 

dx 


(I 


-  S, 


(3) 


Now  imagine  that  (here  arc  two  variables  instead  of  one  For  instance  for  the  tokamak. 
say  that  both  density  and  temperature  gradients  drive  instability  and  that  there  is  both 
anomalous  diffusion  and  anomalous  thermal  conduction  Then  the  equations  which  describe 
the  system  are  heat  and  particle  balance 
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o-  f 

dx  dx 


♦  S„ 


marginal  stability 


and  some  relation  between  I)  and  K 


dT 

dx 


♦  a 


dn 

dx 


P 


(4) 

(5) 

(61 


D  -  {K. 

which  might  come  from  for  instance  Quasi-lincar  theory  Also  we  assume  {  >  0 
lions  can  be  solved  for  K  as  follows  Integrate  Fq  (4)  from  zero  to  r  assuming 
r  -  o  to  get 


(7) 


These  cqua- 
dT 


dx 


-  0  at 


Then  solve  for  D  and 
The  result  is 


fl-Js, 

dx  Jn 


dx 


(8) 


dn 

dx 


from  Eqs  (6*  and  (7).  insert  in  Eq  (5)  and  integrate  from  0  to  r 


K--^J's,d*  ^£sndx 

Then  once  Eq  (9)  is  solved  for  K.  we  can  solve  Eq  (4)  for  T  in  the  standard  way 


(9) 


Let  us  now  examine  the  conditions  under  which  a  solution  to  Eqs  (4-7)  can  be  formed 
In  the  normal  tokamak  configuration  where  the  electrons  are  in  the  plateau  or  banana  regime. 
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negative  —  and  -  both  drive  instability  Therefore  a  and  &  have  the  same  sign  and  boih 

<J\  0\ 

have  opposite  m gn  to  P  Since.  S, .  S„  and  y  are  all  positive,  K  is  greater  lhan  zero  so  that  an  a 
solution  exists  Nov*  imagine  lhai  say  ihc  temperature  gradient  is  stabilizing  ralher  than  desta¬ 
bilizing.  tbai  is  ii  and  P  now  have  the  same  sign  Then,  in  order  to  have  K  >  0. 


. i J- 

f*  S.,dx|  > 

D  f  SIdX 

Hv 

«*o 

p  Jit  1 

(10) 

It  inequality  (101  is  violated  in  some  regions  of  x,  the  plasma  will  be  stable  and  a  marginally 
stable  solution  docs  not  exist  there  In  these  regions  classical  transport  will  apply 


Of  course  ralher  than  solving  equations  like  (5)-<7).  a  more  practical  procedure  is  usually 
to  use  a  numerical  transport  code  in  which  transport  coefficients  jump  to  some  large  value  at 
the  stability  threshold,  and  in  which  the  proper  relation  between  the  different  transport 
coefficients  is  also  used  \x  long  as  these  tumps  in  transport  coefficients  arc  handled  in  a  way 
which  is  numerically  stable,  the  o*dc  will  automatically  handle  the  transition  not  only  between 
stable  and  unstable  plasma  but  also  between  different  instabilities 


For  the  next  stage  of  the  hcirarchy  one  must  know  the  nonlinear  limit  to  the  fluctuation 
level  However  let  us  assume  that  once  the  fluctuation  level  is  known,  the  transport 
coefficients  can  all  be  calculated,  for  instance  via  quasi-lincar  theory  Actually  this  is  not  as  un¬ 
reasonable  as  it  sounds  Sav  that  the  instability  is  driven  by  electrons  and  is  stabilized  by  some 
nonlinear  effect  on  ions  Then  quasi-linear  transport  coefficients  for  the  electrons  should  be 
valid  If  this  be  the  ^ase.  at  least  some  of  the  ton  transport  coefficients  can  usually  be  calculat¬ 
ed  by  invoking  global  conservation  relations  (for  instance  energy  or  momentum  conservation, 
pressure  balance,  etc)  Thus  ihc  third  level  of  our  hcirarchy  assumes  that  transport  coefficients 
can  all  be  calculated  in  terms  of  a  fluctuation  level  which  is  specified  by  invoking  non-linear 
theory,  and/or  experiment,  and/or  numerical  simulation  of  the  instability.  There  arc  at  least 
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iwo  reasons,  why  the  lluctuaiion  level  might  be  required  First,  the  system  may  he  so  strongly 
driven  that  nonlinear  effects  limit  the  fluctuations  to  a  value  smaller  than  that  required  to  main¬ 
tain  a  marginally  stable  profile  Second,  there  may  be  so  many  anomalous  transport  cfTccts  that 
there  is  no  simple  external  mechanism  driving  the  plasma  toward  instability  and  no  simple  re¬ 
laxation  to  a  stable  state  because  of  the  instability 

This  is  the  approach  recently  used  at  the  Naval  Research  Laboratory  to  describe 
anomalous  absorption  and  flux  limitation  in  a  laser  produced  plasma  *  F.arlicr  calculations  like 
this  also  examined  hydrodynamic  flow  in  an  ionosphere  which  has  been  violently  perturbed."1 
and  the  implosion  and  post  implosion  phase  of  theta  pinches  "  Since  these  results  all  require 
a  fluctuation  level  which  is  the  result  of  a  nonlinear  calculation  they  arc  less  reliable  than 
results  obtained  via  calculations  on  the  two  lower  levels  of  the  hcirarchy  which  we  have  been 
describing 

The  fourth  level  of  the  hcirarchy,  and  the  last  one  which  we  will  describe,  is  that  where 
both  the  nonlinear  fluctuation  level  and  the  transport  coefficient  arc  the  results  of  nonlinear  cal¬ 
culations  Needless  to  see  results  on  this  level  arc  still  more  speculative  than  those  on  the 
three  lower  levels  We  will  now  describe  three  calculations  done  at  the  Naval  Research  Labora¬ 
tory  which  illustrate  the  three  levels  They  arc  calculations  of  electron  and  ion  temperature 
profiles  in  tokamaks4.  calculation  of  the  structure  of  transverse  resistive  shocks',  and  calcula¬ 
tions  of  anomalous  absorption  and  thermal  energy  flux  limitation  in  a  laser  produced  plasma  " 

In  a  lokamak.  the  electron  temperature  is  limited  by  some  anomalous  process,  generally 
thought  to  be  instabilities  We  have  attempted  to  study  this  with  the  use  of  a  simple  transport 
code  *  Specifically  our  code  solves  the  equations 

-  f  s?""1"  * K-’ir +  E  •  - >'■  5 - (T- - T>  »" 
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3  1  d  dT,  m  „  T  . 

— n  — - —  rk  — —  +  3vf  —  n  (Tc  -  T, )  (12) 

2  dt  r  dr  dr  M 

E  -  J/«r  03) 

where  Te .  arc  the  electron  (ion)  temperature,  i  is  the  current  density,  E  is  the  toroidal  electric 

field,  assumed  to  he  independent  of  radius.  <r  is  the  electrical  conductivity. 


The  current  density  is  assumed  to  always  he  at  its  steady  slate  value,  so  magnetic  diffusion 
is  neglected  Eq  1 1 1  > - <  1 3 >  are  solved  subject  to  boundary  conditions  that  T,  —  T,  —  lOev  at 
the  limiter,  and  that  the  total  lurrent  is  specified  Every  term  in  Eq  <  1 1 ) - <  1 3 )  is  given  by 
neoclassical  theory  escept  for  one.  the  anomalous  electron  thermal  conduction  Kjn  If  the  plas¬ 
ma  is  unstable.  K,„  is  taken  as  the  Rohm  diffusion  coefficient  K„  In  our  study  the  plasma  is 

rB, 

subject  to  one  of  two  instabilities,  the  internal  kink  tearing  mode  whenever  q  —  — —  "  less 

RBi 

than  unity,  and  the  universal  drill  wave-trapped  electron  mode  whenever  the  shear  strength  H 

is  below  a  critical  value  «  ,  I  he  value  of  «  ,  comes  from  the  linear  theory  of  the  mode  We 
will  not  elaborate  here,  but  a  full  discussion  can  be  found  in  Ref  4  Thus  the  functional  form 

for  k,„  is 


k  *  ”  k  n 


I  ♦  q 


(H'H.t  " 

I  ♦<«/»  ,) 


(14) 


where  n  and  n  arc  targe  integers,  so  that  k,„  turns  on  abruptly  as  the  plasma  becomes  un¬ 
stable  Since  the  quantities  in  the  brackets  change  abruptly  from  zero  to  one.  we  have  called 
them  thermostat  functions  Let  us  now  re-emphasize  that  the  nonlinear  theory  of  the  instabili¬ 
ty  plays  no  role  in  the  analvsis  The  only  thing  needed  to  solve  Fqs  (II)-(I3)  is  the  instability 
condition  which  comes  onlv  from  bmuir  theory 


Now.  knowing  the  profile,  one  can  utilize  quasi-linear  theory  to  calculate  both  the 
anomalous  thermal  conduction,  fluctuation  wave  number  and  fluctuation  amplitude  However, 
even  if  the  behavior  is  governed  by  nonlinear  theory,  the  only  thing  which  will  be  wrong  is  the 
calculation  of  fluctuation  level,  not  the  profile 
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We  Mill  no*  show  some  results  of  (he  calculation  Since  (here  was  nol  much  available 
data  Irons  PLT  a(  (he  time,  see  did  onl>  one  run  The  parameter*  are  H  -  35  KCi.  /c,  -  4.  a 
central  density  of  n0  “  4  2*10"  and  a  hydrogen  plasma  The  electron  and  ion  temperature,  q. 
and  thermostat  functions  for  q  and  M  are  given  in  Fig  3  Clearly  the  plasma  is  in  a  marginally 

stable  state  for  internal  kink -tearing  modes  lor  0<  -<0  4,  and  to  trapped  electron  modes  for 

a 

0  4  <  -  <0  f  The  outer  region  is  classical 
a 

Another  tokamak  which  we  have  simulated  is  TFR  In  TFR  more  than  half  of  the  input 
power  is  generally  radiated  away  by  impurity  radiation  However,  as  is  clearly  indicated  in 
Ref  13.  most  of  this  is  oxygen  line  radiation  arising  from  the  plasma  edge  where  the  tempera¬ 
ture  is  low  enough  that  the  plasma  is  stable  Thus  the  question  is  how  does  the  energy  get 
from  the  center  to  the  edge  Reference  13  also  shows  that  destruction  of  magnetic  surfaces  due 

to  internal  kink  and  tearing  modes  is  insufficient  to  account  for  this  energy  transport  Here  we 
examine  whether  drift  and  trapped  particle  instabilities  can  provide  the  remainder  of  the  energy 
transport 

One  very  interesting  experimental  result  for  TFR"  is  that  the  temperature  half  width 
2ir  IT  >  depends  on  only  a  single  parameter,  q  (a  >.  even  though  two  parameters.  B  and  I  arc  in¬ 
dependently  varied  Our  calculation  basically  confirms  this  result  In  Fig  4  is  shown  (solid 
line*  the  predictions  by  our  code  for  AHT  >  as  a  function  of  q  ‘(a  • 

Here  we  have  assumed  7r,  -  3  and  n„  -  b  *10"  Various  dots  shown  the  calculated 
points  Although  there  is  some  scatter  to  the  points,  Ar  <T  i  basically  does  depend  on  the  single 
parameter  q  (a  I  The  dashed  line  shows  the  experimental  result  from  TFR 

Another  interesting  result  from  the  TFR  experiment  is  that  the  position  of  the  q  -  I  sur¬ 
face.  as  defined  by  the  radial  position  of  the  node  in  the  saw  tooth  oscillations  in  the  soft  x-ray 
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signal,  also  depends  only  on  the  single  parameter  4  (a  I  In  I  ig  4  the  dash  dot  line  shows  that 
the  predicted  value  for  the  <4  -  1  surface,  as  defined  by  the  radial  position  where  the  4  ther¬ 
mostat  function  drops  to  halt  its  maximum  value.  also  depends  principally  on  this  single  param¬ 
eter 

In  I  ig  >(a»  and  (hi  are  predictions  of  the  code  lor  electron  temperature  compared  to  ex¬ 
periment,  for  three  discharges  having  I  -  140K.A,  and  B  -  25,  40  and  50  Kg  Notice  that  the 
code  does  give  good  agreement  with  experiment  as  regards  temperature  profile  In  |-'ig  6  arc 
plotted  the  radial  dependence  of  T  1  .  4.  for  the  run  with  B  -  50  Kg.  1  —  140  KA 

Other  4uantities  ol  interest  arc  the  k.  ol  the  fluctuation  and  predicted  values  of  c  i  T,  as 
a  function  of  radius  Ihese  ate  shown  in  I  ig  7  for  a  discharge  with  60  KV.  300  KA  The  fluc¬ 
tuation  level  seems  to  he  comparable  to  that  recently  reported  4  e  </>  !f  =  (3-41  x  |0 

In  addition  to  calculating  the  temperature  profile  for  TI  R  numerous  calculations  wetc 

also  done  for  Alcator  Perhaps  the  most  striking  experimental  result  fot  Alcator  is  the  increase 

of  energv  confinement  time  with  density  Our  theory  does  predict  this  haste  dependence  In 

I  ig  H  is  shown  the  calculated  dependence  of  energy  confinement  time  on  central  density  fot 

three  choices  of  field  and  current  B  —  Ml  k(i.  I  —  1  (HI  kA.  B  —  75  kCi.  I  “  100  k.A,  and  B 

-  75  KG.  I  -  150  k  A  In  all  cases.  /  —  I  and  we  used  a  hydrogen  plasma  Notice  that  for 

central  densities  less  than  about  2  «  10  4  cm  .  the  confinement  time  increases  roughly  linearly 

with  densitv  I  or  larger  n  the  confinement  time  begins  to  decrease  again 

When  the  densitv  is  suftiucntlv  high  however,  the  plasma  is  in  the  Pfirsch  Schlutcr  re¬ 
gime.  and  there  is  no  trapped  particle  instability  Then  the  energy  confinemeni  lime  begins  10 

decrease  with  densitv  lor  all  three  choices  of  current  and  field  the  points  with  n0  —  5  x  |0U 
were  completely  stable  Not  only  arc  there  no  trapped  particle  instabilities  in  the  region  of 
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maximum  gradient,  q  was  everywhere  larger  than  one  so  that  arc  no  MHO  modes  m  the  center 

cither  Apparently  the  largest  confinement  time  occurs  just  on  the  transition  from  anomalous 

to  classical  thermal  conduction  The  highest  electron  temperature  however  are  in  the  regime  of 

anomalous  thermal  conduction  and  low  confinement  time  In  f  ig  *Ja  and  b  are  plotted  the 
radial  dependence  of  Tc  and  T  for  different  densities  for  the  case  B  -  75  k(i.  I  -  100  k\ 

Notice  that  at  low  density,  where  the  thermal  conduction  is  anomalous,  the  temperature  profiles 
are  quite  peaked  At  higher  density,  where  the  thermal  conduction  is  classical,  the  temperature 
profiles  are  quite  broad 

To  summarize,  our  calculations  of  tokamak  temperature  profiles  do  give  qualitative  and 

even  some  quantitative  agreement  with  measured  tokamak  temperature  profiles  The  theory  is 

simple  in  that  nonlinear  theory  of  the  instability  never  enters  in  However  there  arc  no  simple 

scaling  laws  which  explain  tokamak  behavior,  instead  the  temperature  profiles  result  only  from 
a  numerical  solution  of  the  electron  and  ion  energy  equations  with  the  appropriate  anomalous 

electron  thermal  conductivity 

Now  we  consider  a  different  plasma  configuration,  a  cross  field  collisionlcss  resistive 

shock  Here,  as  we  will  see.  marginal  stability  appears  to  be  a  viable  concept,  but  two  transport 
processes,  electron  <zrr</  ion  heating  plav  important  roles  In  a  shock,  fluid  convection  (that  is 

the  \  term  in  the  momentum  equation!  tends  to  steepen  the  dcnsitv  profile  Since  the 
dx 

magnetic  field  is  frozen  into  the  flow  it  also  steepens  and  the  current  thereby  increases  How¬ 
ever  at  some  point,  the  current  becomes  so  great  that  ion  acoustic  waves  arc  driven  unstable 
These  waves  then  grow,  the  anomalous  resistivity  increases,  and  the  shock  profile  will  then 
tend  to  broaden  The  shock  profile  will  then  be  determined  by  the  condition  that  ton  acoustic 
waves  arc  everywhere  marginally  stable,  that  is' 


c  dB 

Tt 

I  : 

1  a. 

1  wfjll 

i  : 

4nnc  dx 

1  ♦ 

1  12  T, 
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From  ihc  Rankine-Hugimst  telaiions,  ihc  magnetic  held  density  and  temperature  are  known 


are  all  known  upstream  and  downstream  Thus,  once  »t  know  — ,  the  ion  to  electron  tern- 

•  C 


perature  ratio.  »c  van  estimate  IB  <upsircam)-B  (downstreaml|/L,  and  thereby  estimate 

ds 

the  shock  width  L, 


The  problem  now  i\  to  estimate  the  ion  to  electron  temperature  ratio  Tor  an  ion  acoustic 
wase  at  marginal  stability,  this  van  be  done  bv  vomparing  the  heating  rates  of  ions  and  elec¬ 
trons  Say  that  the  ion  wase  transfers  momentum  from  tons  to  electrons  at  a  rate  P  Then  it 
can  be  shown,  from  resonant  quasi-lincar  theory,  that  in  the  reference  frame  in  which  the  ions 
arc  at  rest,  enctgv  is  transferred  from  electrons  to  ions  at  a  rate  (q»/k  IP.  where  lu  k  I  is  Ihc 
phase  velocity  of  the  ion  acoustic  wave  If  the  electron  drift  velocity  is  denoted  u.  then 
momentum  and  energy  conservation  equations  for  electrons  and  ions  read 


nmu  -  P 

(a  1 

(16) 

nmuu  *  1  ?  n  T(  - 

“  p  _ 

k 

JnT 

* 

(b) 

I  rom  Tqs  t|f>a  and  hi.  it  is  a  simple  matter  to  solve  for  the  ratio  of  heating  rates. 


1 

T 


I  u  <U  k  • 
tid  k  t 


(171 


If  both  electrons  and  ions  arc  substantial^  healed  bv  the  shock,  the  temperature  ratio  should  be 
roughly  equal  to  the  ratio  of  heating  rates  or 


U  u  k  I 

w  k  v\2 


I  *• 


i  ♦  i: 


t, 

T, 


l*M  i 


c*p 


T, 

T 


1*1*12 


T 


t  • 


T, 

-  I  =  — -<  1 8 1 


T,  T, 

I  quation  (IRi  is  a  tranccniental  equation  for  Tor  a  hydrogen  plasma,  we  find  — -  —  7.5 

*  •  t 

T 

Then  knowing  .  Iq  tl't  gives  a  simple  estimate  for  shock  width  Tor 

"  i 
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—  =  70.  — =  7  5-/^"  0  01.  (19) 

<•<.,  4irnc  L,  a/  M  BJ 

and  a  Mach  two  shock,  this  gives  the  result  L,  — 10c  which  is  consistent  with  the  transverse 
shock  experiments  done  at  Culham  laboratory  From  this  estimate  of  shock  width,  one  can 
calculate  a  resistivity  and  from  the  resistivity,  a  fluctuation  level  This  works  out  to  he  about 


—  5*10  1  (20) 

C  <t> 

let  us  rc  cmphasi/c  that  at  no  point  in  Thi*  calculation  an  estimate  of  -  from  nonlinear 

I  c 

theory  ever  required 


\ctually.  one  can  do  much  better  than  estimate  these  quantities  In  Ref  3.  fluid  equa¬ 
tions  for  n.  Tf.  T  and  B.  coupled  to  wave  equations  for  lb  (or  64)  ion  acoustic  fluctuations  r at 
different  wave  vcviorl  arc  numericallv  integrated  from  upstream  to  downstream  The  waves  arc 
assumed  to  grow  or  damp  at  the  local  linear  growth  or  damping  rate,  and  quasi-lincar  theory  is 
used  to  calculate  resistivity  as  well  as  electron  and  ion  heating  rate  It  was  found  that  the  shock 

profile  did  in  fact  remain  at  a  marginally  stable  current  for  nearly  the  entire  shock  profile 
Results  of  such  a  calculation  arc  shown  in  Fig  10  where  spatial  profiles  of  n  B  Tr.  T  and 

are  shown  To  summarize,  the  structure  of  transverse  resistive  shocks  in  hydrogen 

seems  to  be  consistent  with  ion  acoustic  waves  being  at  marginal  stability  everywhere  in  the 
profile  I  stimates  and  calculations  of  shock  width  and  fluctuation  level  arc  in  good  quantitative 
agreement  with  experiments 


As  a  final  example,  we  consider  the  absorption  of  laser  light  in  a  laser  produced  plasma 
This  is  discussed  much  more  fullv  in  Ref  9  and  the  description  here  will  be  very  brief  The 
idea  is  to  use  a  numerical  solution  of  the  fluid  and  wave  equations  with  anomalous  transport 
There  are  two  instabilities  which  provide  this  anomalous  transport  First,  there  is  Brillouin 
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backscatter .  where  in  the  laser  light  decays  into  a  reflected  wave  and  an  ion  acoustic  wave,  and 
second,  there  is  an  ion  acoustic  instability  generated  b>  the  return  current 

For  Bnllouin  backscatter.  if  the  growth  rate  in  homogeneous  media  is  y0.  then  the  spatial 
amplification  in  inhomogeneous  for  cither  an  undamped.'  or  strongly  damped  ion  acoustic 
wave,  is 

s  :>  2 

•  yoL  pi, 

C\? 

where  l  is  the  si/c  of  the  region  of  phase  coherence  between  incident,  reflected,  and  ion 
acoustic  wave  1  hi\  amplification  turns  out  to  be  so  large  that  any  theory  based  upon  it  will 
gise  nearly  total  backsc.iticr  in  vntuallv  am  circumstance  However  when  examining  the  ampli¬ 
tude  of  the  ion  acoustic  wave  gcncraied  one  secs  that  it  is  very  large,  so  large  that  any  ion 
trapping  should  be  verv  important  and  linear  theory  of  the  instability  invalid  Thus  a  descrip¬ 
tion  of  the  process  has  to  be  based  on  a  nonlinear  description  of  the  instability  1  arlicr  theory 
"  :  has  shown  that  if  trapping  is  important  one  wav  the  instability  can  be  modeled  is  by  reduc¬ 
ing  the  growth  rate  >  b»  a  factor  of  between  about  five  and  ten  We  have  adopted  such  an  ap¬ 
proach  to  model  the  effect  of  Hnlloum  back  vat  ter  Thus  we  need  not  only  a  nonlinear  theory 
of  the  fluctuation  amplitude  but  also  a  nonlinear  reduction  in  the  reflection  due  to  stimulated 
Brilloutn  backveatter 

The  second  mstahilitv  which  we  consider  is  the  ion  acoustic  instability  driven  by  a  return 

current  If  electrons  conduct  heat  but  cam  no  current,  a  flux  of  energetic  particles  in  one 
direction  is  balanced  bv  a  flux  of  low  velocitv  particles  going  the  other  way  This  return  current 

can  drive  ion  acoustic  waves  unstable  There  arc  three  principle  effects  of  this  instability,  first 

the  electron  thermal  conductivity  is  reduced,  second,  there  is  an  electron  ion  energy  exchange. 

and  third,  there  is  anomalous  absorption  of  laser  light  resulting  from  the  scattering  of  laser  light 
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on  the  ion  density  fluctuations  A  linear  and  guasi- linear  theory  of  this  instability  has  been 
published  recently,*1  and  a  self  consistent  steady  state  treatment  of  the  coupling  of  anomalous 
absorption  and  flux  limitation  has  also  been  presented  "  Mere,  a  fluctuation  lescl  is  required, 
but  then  with  this  fluctuation  level,  the  transport  coefficients  are  related  by  (Juasi-lincar  theory 

The  approach  in  Ref  is  to  numerically  solve  the  fluid  equations  where  anomalous  tran¬ 
sport  and  backscatter  results  from  ihcsc  two  instabilities  \k  c  have  simulated  three  types  of  Nd 
laser  pulses  hirst  single  short  l7t)p  seel  pulses  on  a  target  which  has  an  initial  density  gradient 
scale  length  of  I0»i  This  simulates  many  of  the  short  pulse  experiments,  especially  on  slab  tar¬ 
gets  We  find  that  the  absorption  is  principally  resonant  absorption,  but  it  gets  a  strong  boost 
from  return  current  driven  turbulence  from  irradianccs  of  10  4  \k  cm'  to  .1*10"'  W  /cm  ;  the 
absorption  is  about  40*.  with  about  2Vk  resonant  absorption  lalso  shown  in  I  ig  II  is  back- 
waiter  and  specular  reflection  vs  irradiancc  I  This  is  in  reasonable  agreement  with  data  taken 
from  mans  laboratories  In  f  ig  12  arc  shown  expansion  velocity,  electron  temperature  and  flux 
limit.  Q  nm\  as  a  function  of  irradiancc  The  points  are  measured  electron  temperatures 
from  the  NRl  experiment 

A  second  simulation  is  of  a  double  structured  pulse  Recent  experiments  at  NRL  have 
shown  that  if  the  main  laser  pulse  illuminates  a  prcpulsc  formed  plasma,  vers  strong  backscatter 
results  *  Ike  model  this  bv  assuming  a  100»im  scale  length  for  n  <0  In,,  and  a  10>j  scale 
length  for  N  >\,  and  illuminate  with  a  single  70p  sec  pulse  Shown  in  f  ig  13  is  absorption. 

backscatter.  resonant  absorption  and  specular  'cflection  as  a  function  of  irradiancc  Also  shown 
arc  experimental  measurements  of  backscatter  as  a  function  of  irradiancc  Clearly  there  is  good 

qualitative  agreement  betwen  theory  and  experiment  here 

Finally  we  have  simulated  long  pulse  experiments  also  Mere  we  assume  a  gradient  scale 
length  of  IOOjj  and  take  a  pulse  which  rises  up  to  a  final  steady  irradicncc  and  then  remains 
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constant  In  I  igs  14  and  IS  arc  shown  plots  til  absorption,  back  scat  ter  and  specular  reflection 

versus  irradiance.  and  then  plots  ol  clcvtron  temperature,  expansion  velocity  and  flue  limit 

versus  irradiance  1  he  points  on  the  figures  arc  experimental  measures  of  absorption  ■'  •’*  for 

long  pulse  experiments  To  summarize,  the  fluid  simulation  with  anomalous  transport  seems  to 

provide  qualitative  and  even  lair  quantitative  agreement  with  experiments  on  laser  produced 
plasmas  fin  a  vanciv  of  laser  pulse  shapes  and  over  nearl>  five  orders  of  magnitude  in  irradi- 

ance 


Now  what  about  our  initial  question  concerning  whether  simple  scaling  laws  can  be 
derived  from  first  principles’  The  entire  thrust  of  this  paper  is  that  they  really  cannot  Howev¬ 
er  this  does  not  mean  that  no  progress  can  he  made,  on  the  contrary  it  seems  that  one  can  go 
tar  with  fluid  simulations  with  the  instahilitv  entering  via  anomalous  transport  In  some  cases. 

all  that  is  needed  for  good  modeling  is  the  linear  stahilitv  threshold  In  other  cases  nonlinear 
thcorv  and/or  relation  between  different  transport  coefficients  is  needed  However  in  all  cases. 

the  scaling  comes  from  a  numerical  solution  of  the  fluid  equations,  not  from  a  simple  scaling 

law 
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INTERFACE  WITH  FLUID  CODES 


I  SOME  PARAMETER 

FOR  A  WIDE  RANGE  OF  EXTERNAL  SOURCE 
STRENGTHS,  FLUID  CODE  WILL  AUTOMATICALLY  PICK 
OUT  A  MARGINALLY  STABLE  PROFILE. 


I  tg|  I  —  \  pnwhlc  p*<*«  i»f  transport  inefficient  vercuv  w«mc  relevant  parameter 


r 


MU  Ml  MUKAMH  M  Rfr PORT  4066 


TURBULENCE 

TRANSPORT 

LEVEL 

COEFFICIENT 

PROFILE 


CONVENTIONAL  NONLINEAR  THEORY 
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q(a)"1 


T»f  4  —  Temperature  half  •  nit h  and  •  •dth  of  the  ffr/  •  I  region  in  TUI  a*  a  function  of  $t<il  1  The  aolid 
and  dath  dot  line  are  apcrnumaie  At*  to  the  numernalU  obtained  data  point*  The  da*hed  line  i%  the  expert 
mental!*  observed  temperature  half  •  nlth 
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SPECTRUM 


F.ncrgy  confinement  lime  *e«u*  centril  dentil*  in  Akatnr 


ft*  «  —  Calculated  electron  <•  >  end  ton  Ibi  temperature  prnAlet  at 
different  denaitiet  Man  a  ploi  tel  of  meaaured  ce*''*?'  electron  and 
ton  temperature*  *t  denattt 
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Nd  SHORT  PULSE 
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(  aWuUtcit  U«cr  Itfht  tbvifpiion  rc%**n*m  alnorplion  hdckvjHff  ind  \pct ulir  refWciion  vertut  irrsdisnce 
for  a  \d  I ttsct  plasma  •  iih  a  pultc  of  **0  pv»  f«hm 
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Nd  DOUBLE-STRUCTURED  (30-1<XM 
SHORT  PULSE 
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13  —  C  akuUtcd  iHvKpOon.  baik\*.«tter  remnant  absorption  and  specular  reflection 
for  a  double  structure  Nd  laser  pulse 
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